Monozygotic twins discordant for schizophrenia differ in maturation and synaptic transmission.

Schizophrenia affects approximately 1% of the world population. Genetics, epigenetics, and environmental factors are known to play a role in this psychiatric disorder. While there is a high concordance in monozygotic twins, about half of twin pairs are discordant for schizophrenia. To address the question of how and when concordance in monozygotic twins occur, we have obtained fibroblasts from two pairs of schizophrenia discordant twins (one sibling with schizophrenia while the second one is unaffected by schizophrenia) and three pairs of healthy twins (both of the siblings are healthy). We have prepared iPSC models for these 3 groups of patients with schizophrenia, unaffected co-twins, and the healthy twins. When the study started the co-twins were considered healthy and unaffected but both the co-twins were later diagnosed with a depressive disorder. The reprogrammed iPSCs were differentiated into hippocampal neurons to measure the neurophysiological abnormalities in the patients. We found that the neurons derived from the schizophrenia patients were less arborized, were hypoexcitable with immature spike features, and exhibited a significant reduction in synaptic activity with dysregulation in synapse-related genes. Interestingly, the neurons derived from the co-twin siblings who did not have schizophrenia formed another distinct group that was different from the neurons in the group of the affected twin siblings but also different from the neurons in the group of the control twins. Importantly, their synaptic activity was not affected. Our measurements that were obtained from schizophrenia patients and their monozygotic twin and compared also to control healthy twins point to hippocampal synaptic deficits as a central mechanism in schizophrenia.

Activity-Dependent Transcriptional Program in NGN2+ Neurons Enriched for Genetic Risk for Brain-Related Disorders.

BACKGROUND: Converging evidence from large-scale genetic and postmortem studies highlights the role of aberrant neurotransmission and genetic regulation in brain-related disorders. However, identifying neuronal activity-regulated transcriptional programs in the human brain and understanding how changes contribute to disease remain challenging. METHODS: To better understand how the activity-dependent regulome contributes to risk for brain-related disorders, we profiled the transcriptomic and epigenomic changes following neuronal depolarization in human induced pluripotent stem cell-derived glutamatergic neurons (NGN2) from 6 patients with schizophrenia and 5 control participants. RESULTS: Multiomic data integration associated global patterns of chromatin accessibility with gene expression and identified enhancer-promoter interactions in glutamatergic neurons. Within 1 hour of potassium chloride-induced depolarization, independent of diagnosis, glutamatergic neurons displayed substantial activity-dependent changes in the expression of genes regulating synaptic function. Depolarization-induced changes in the regulome revealed significant heritability enrichment for schizophrenia and Parkinson's disease, adding to mounting evidence that sequence variation within activation-dependent regulatory elements contributes to the genetic risk for brain-related disorders. Gene coexpression network analysis elucidated interactions among activity-dependent and disease-associated genes and pointed to a key driver (NAV3) that interacted with multiple genes involved in axon guidance. CONCLUSIONS: Overall, we demonstrated that deciphering the activity-dependent regulome in glutamatergic neurons reveals novel targets for advanced diagnosis and therapy.

Precise Therapeutic Targeting of Distinct NRXN1+/- Mutations.

As genetic studies continue to identify risk loci that are significantly associated with risk for neuropsychiatric disease, a critical unanswered question is the extent to which diverse mutations--sometimes impacting the same gene-- will require tailored therapeutic strategies. Here we consider this in the context of rare neuropsychiatric disorder-associated copy number variants (2p16.3) resulting in heterozygous deletions in NRXN1, a pre-synaptic cell adhesion protein that serves as a critical synaptic organizer in the brain. Complex patterns of NRXN1 alternative splicing are fundamental to establishing diverse neurocircuitry, vary between the cell types of the brain, and are differentially impacted by unique (non-recurrent) deletions. We contrast the cell-type-specific impact of patient-specific mutations in NRXN1 using human induced pluripotent stem cells, finding that perturbations in NRXN1 splicing result in divergent cell-type-specific synaptic outcomes. Via distinct loss-of-function (LOF) and gain-of-function (GOF) mechanisms, NRXN1+/- deletions cause decreased synaptic activity in glutamatergic neurons, yet increased synaptic activity in GABAergic neurons. Stratification of patients by LOF and GOF mechanisms will facilitate individualized restoration of NRXN1 isoform repertoires; towards this, antisense oligonucleotides knockdown mutant isoform expression and alters synaptic transcriptional signatures, while treatment with ß-estradiol rescues synaptic function in glutamatergic neurons. Given the increasing number of mutations predicted to engender both LOF and GOF mechanisms in brain disease, our findings add nuance to future considerations of precision medicine.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome.

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

Multi-omic profiling of the developing human cerebral cortex at the single-cell level.

The cellular complexity of the human brain is established via dynamic changes in gene expression throughout development that is mediated, in part, by the spatiotemporal activity of cis-regulatory elements (CREs). We simultaneously profiled gene expression and chromatin accessibility in 45,549 cortical nuclei across six broad developmental time points from fetus to adult. We identified cell type-specific domains in which chromatin accessibility is highly correlated with gene expression. Differentiation pseudotime trajectory analysis indicates that chromatin accessibility at CREs precedes transcription and that dynamic changes in chromatin structure play a critical role in neuronal lineage commitment. In addition, we mapped cell type-specific and temporally specific genetic loci implicated in neuropsychiatric traits, including schizophrenia and bipolar disorder. Together, our results describe the complex regulation of cell composition at critical stages in lineage determination and shed light on the impact of spatiotemporal alterations in gene expression on neuropsychiatric disease.

Better together: Non-additive interactions between schizophrenia risk genes.

Zhang, Zhang, Forrest et al combine allele-specific open chromatin (ASoC) mapping and CRISPR-editing to evaluate the functional impact of schizophrenia risk variants on human neuronal gene expression, synaptic development, and function. In doing so, they uncover surprising non-additive effects between target genes regulated by the same risk variant.

Schizophrenia-associated somatic copy-number variants from 12,834 cases reveal recurrent NRXN1 and ABCB11 disruptions.

While germline copy-number variants (CNVs) contribute to schizophrenia (SCZ) risk, the contribution of somatic CNVs (sCNVs)-present in some but not all cells-remains unknown. We identified sCNVs using blood-derived genotype arrays from 12,834 SCZ cases and 11,648 controls, filtering sCNVs at loci recurrently mutated in clonal blood disorders. Likely early-developmental sCNVs were more common in cases (0.91%) than controls (0.51%, p = 2.68e-4), with recurrent somatic deletions of exons 1-5 of the NRXN1 gene in five SCZ cases. Hi-C maps revealed ectopic, allele-specific loops forming between a potential cryptic promoter and non-coding cis-regulatory elements upon 5' deletions in NRXN1. We also observed recurrent intragenic deletions of ABCB11, encoding a transporter implicated in anti-psychotic response, in five treatment-resistant SCZ cases and showed that ABCB11 is specifically enriched in neurons forming mesocortical and mesolimbic dopaminergic projections. Our results indicate potential roles of sCNVs in SCZ risk.

Schizophrenia Risk Mapping and Functional Engineering of the 3D Genome in Three Neuronal Subtypes.

Common variants associated with schizophrenia are concentrated in non-coding regulatory sequences, but their precise target genes are context-dependent and impacted by cell-type-specific three-dimensional spatial chromatin organization. Here, we map long-range chromosomal conformations in isogenic human dopaminergic, GABAergic, and glutamatergic neurons to track developmentally programmed shifts in the regulatory activity of schizophrenia risk loci. Massive repressive compartmentalization, concomitant with the emergence of hundreds of neuron-specific multi-valent chromatin architectural stripes, occurs during neuronal differentiation, with genes interconnected to genetic risk loci through these long-range chromatin structures differing in their biological roles from genes more proximal to sequences conferring heritable risk. Chemically induced CRISPR-guided chromosomal loop-engineering for the proximal risk gene SNAP91 and distal risk gene BHLHE22 profoundly alters synaptic development and functional activity. Our findings highlight the large-scale cell-type-specific reorganization of chromosomal conformations at schizophrenia risk loci during neurodevelopment and establish a causal link between risk-associated gene-regulatory loop structures and neuronal function.

R-loop landscapes in the developing human brain are linked to neural differentiation and cell-type specific transcription.

Here, we construct genome-scale maps for R-loops, three-stranded nucleic acid structures comprised of a DNA/RNA hybrid and a displaced single strand of DNA, in the proliferative and differentiated zones of the human prenatal brain. We show that R-loops are abundant in the progenitor-rich germinal matrix, with preferential formation at promoters slated for upregulated expression at later stages of differentiation, including numerous neurodevelopmental risk genes. RNase H1-mediated contraction of the genomic R-loop space in neural progenitors shifted differentiation toward the neuronal lineage and was associated with transcriptomic alterations and defective functional and structural neuronal connectivity in vivo and in vitro. Therefore, R-loops are important for fine-tuning differentiation-sensitive gene expression programs of neural progenitor cells.

Contributions of circadian clock genes to cell survival in fibroblast models of lithium-responsive bipolar disorder.

Bipolar disorder (BD) is characterized by mood episodes, disrupted circadian rhythms and gray matter reduction in the brain. Lithium is an effective pharmacotherapy for BD, but not all patients respond to treatment. Lithium has neuroprotective properties and beneficial effects on circadian rhythms that may distinguish lithium responders (Li-R) from non-responders (Li-NR). The circadian clock regulates molecular pathways involved in apoptosis and cell survival, but how this overlap impacts BD and/or lithium responsiveness is unknown. In primary fibroblasts from Li-R/Li-NR BD patients and controls, we found patterns of co-expression among circadian clock and cell survival genes that distinguished BD vs. control, and Li-R vs. Li-NR cells. In cellular models of apoptosis using staurosporine (STS), lithium preferentially protected fibroblasts against apoptosis in BD vs. control samples, regardless of Li-R/Li-NR status. When examining the effects of lithium treatment of cells in vitro, caspase activation by lithium correlated with period alteration, but the relationship differed in control, Li-R and Li-NR samples. Knockdown of Per1 and Per3 in mouse fibroblasts altered caspase activity, cell death and circadian rhythms in an opposite manner. In BD cells, genetic variation in PER1 and PER3 predicted sensitivity to apoptosis in a manner consistent with knockdown studies. We conclude that distinct patterns of coordination between circadian clock and cell survival genes in BD may help predict lithium response.

Integrating genetics and transcriptomics to study major depressive disorder: a conceptual framework, bioinformatic approaches, and recent findings.

Major depressive disorder (MDD) is a complex and heterogeneous psychiatric syndrome with genetic and environmental influences. In addition to neuroanatomical and circuit-level disturbances, dysregulation of the brain transcriptome is a key phenotypic signature of MDD. Postmortem brain gene expression data are uniquely valuable resources for identifying this signature and key genomic drivers in human depression; however, the scarcity of brain tissue limits our capacity to observe the dynamic transcriptional landscape of MDD. It is therefore crucial to explore and integrate depression and stress transcriptomic data from numerous, complementary perspectives to construct a richer understanding of the pathophysiology of depression. In this review, we discuss multiple approaches for exploring the brain transcriptome reflecting dynamic stages of MDD: predisposition, onset, and illness. We next highlight bioinformatic approaches for hypothesis-free, genome-wide analyses of genomic and transcriptomic data and their integration. Last, we summarize the findings of recent genetic and transcriptomic studies within this conceptual framework.

The functional and evolutionary impacts of human-specific deletions in conserved elements.

Conserved genomic sequences disrupted in humans may underlie uniquely human phenotypic traits. We identified and characterized 10,032 human-specific conserved deletions (hCONDELs). These short (average 2.56 base pairs) deletions are enriched for human brain functions across genetic, epigenomic, and transcriptomic datasets. Using massively parallel reporter assays in six cell types, we discovered 800 hCONDELs conferring significant differences in regulatory activity, half of which enhance rather than disrupt regulatory function. We highlight several hCONDELs with putative human-specific effects on brain development, including HDAC5, CPEB4, and PPP2CA. Reverting an hCONDEL to the ancestral sequence alters the expression of LOXL2 and developmental genes involved in myelination and synaptic function. Our data provide a rich resource to investigate the evolutionary mechanisms driving new traits in humans and other species.

Convergent coexpression of autism-associated genes suggests some novel risk genes may not be detectable in large-scale genetic studies.

Autism spectrum disorder (ASD) is a heritable neurodevelopmental disorder characterized by deficits in social interactions and communication. Protein-altering variants in many genes have been shown to contribute to ASD; however, understanding the convergence across many genes remains a challenge. We demonstrate that coexpression patterns from 993 human postmortem brains are significantly correlated with the transcriptional consequences of CRISPR perturbations in human neurons. Across 71 ASD risk genes, there was significant tissue-specific convergence implicating synaptic pathways. Tissue-specific convergence was further demonstrated across schizophrenia and atrial fibrillation risk genes. The degree of ASD convergence was significantly correlated with ASD association from rare variation and differential expression in ASD brains. Positively convergent genes showed intolerance to functional mutations and had shorter coding lengths than known risk genes even after removing association with ASD. These results indicate that convergent coexpression can identify potentially novel genes that are unlikely to be discovered by sequencing studies.

Pandemic city: Village-in-a-dish unlocks dynamic genetic effects in the brain.

In this issue, Wells et al. combine genetic multiplexing ("village-in-a-dish") and Stem-cell-derived NGN2-accelerated Progenitors (SNaPs) to evaluate genotype-phenotype relationships across 100 donors in the context of Zika virus infection in the developing brain. This resource will be broadly applicable in uncovering how genetic variation underlies risk for neurodevelopmental disorders.

Non-additive effects of schizophrenia risk genes reflect convergent downstream function.

Genetic studies of schizophrenia (SCZ) reveal a complex polygenic risk architecture comprised of hundreds of risk variants, the majority of which are common in the population at-large and confer only modest increases in disorder risk. Precisely how genetic variants with individually small predicted effects on gene expression combine to yield substantial clinical impacts in aggregate is unclear. Towards this, we previously reported that the combinatorial perturbation of four SCZ risk genes ("eGenes", whose expression is regulated by common variants) resulted in gene expression changes that were not predicted by individual perturbations, being most non-additive among genes associated with synaptic function and SCZ risk. Now, across fifteen SCZ eGenes, we demonstrate that non-additive effects are greatest within groups of functionally similar eGenes. Individual eGene perturbations reveal common downstream transcriptomic effects ("convergence"), while combinatorial eGene perturbations result in changes that are smaller than predicted by summing individual eGene effects ("sub-additive effects"). Unexpectedly, these convergent and sub-additive downstream transcriptomic effects overlap and constitute a large proportion of the genome-wide polygenic risk score, suggesting that functional redundancy of eGenes may be a major mechanism underlying non-additivity. Single eGene perturbations likewise fail to predict the magnitude or directionality of cellular phenotypes resulting from combinatorial perturbations. Overall, our results indicate that polygenic risk cannot be extrapolated from experiments testing one risk gene at a time and must instead be empirically measured. By unravelling the interactions between complex risk variants, it may be possible to improve the clinical utility of polygenic risk scores through more powerful prediction of symptom onset, clinical trajectory, and treatment response, or to identify novel targets for therapeutic intervention.

Stem Cell Models for Context-Specific Modeling in Psychiatric Disorders.

Genome-wide association studies reveal the complex polygenic architecture underlying psychiatric disorder risk, but there is an unmet need to validate causal variants, resolve their target genes(s), and explore their functional impacts on disorder-related mechanisms. Disorder-associated loci regulate transcription of target genes in a cell type- and context-specific manner, which can be measured through expression quantitative trait loci. In this review, we discuss methods and insights from context-specific modeling of genetically and environmentally regulated expression. Human induced pluripotent stem cell-derived cell type and organoid models have uncovered context-specific psychiatric disorder associations by investigating tissue-, cell type-, sex-, age-, and stressor-specific genetic regulation of expression. Techniques such as massively parallel reporter assays and pooled CRISPR (clustered regularly interspaced short palindromic repeats) screens make it possible to functionally fine-map genome-wide association study loci and validate their target genes at scale. Integration of disorder-associated contexts with these patient-specific human induced pluripotent stem cell models makes it possible to uncover gene by environment interactions that mediate disorder risk, which will ultimately improve our ability to diagnose and treat psychiatric disorders.

Induction of dopaminergic neurons for neuronal subtype-specific modeling of psychiatric disease risk.

Dopaminergic neurons are critical to movement, mood, addiction, and stress. Current techniques for generating dopaminergic neurons from human induced pluripotent stem cells (hiPSCs) yield heterogenous cell populations with variable purity and inconsistent reproducibility between donors, hiPSC clones, and experiments. Here, we report the rapid (5 weeks) and efficient (~90%) induction of induced dopaminergic neurons (iDANs) through transient overexpression of lineage-promoting transcription factors combined with stringent selection across five donors. We observe maturation-dependent increase in dopamine synthesis and electrophysiological properties consistent with midbrain dopaminergic neuron identity, such as slow-rising after- hyperpolarization potentials, an action potential duration of ~3 ms, tonic sub-threshold oscillatory activity, and spontaneous burst firing at a frequency of ~1.0-1.75 Hz. Transcriptome analysis reveals robust expression of genes involved in fetal midbrain dopaminergic neuron identity. Specifically expressed genes in iDANs, as well as those from isogenic induced GABAergic and glutamatergic neurons, were enriched in loci conferring heritability for cannabis use disorder, schizophrenia, and bipolar disorder; however, each neuronal subtype demonstrated subtype-specific heritability enrichments in biologically relevant pathways, and iDANs alone were uniquely enriched in autism spectrum disorder risk loci. Therefore, iDANs provide a critical tool for modeling midbrain dopaminergic neuron development and dysfunction in psychiatric disease.

Common genetic variation impacts stress response in the brain.

To explain why individuals exposed to identical stressors experience divergent clinical outcomes, we determine how molecular encoding of stress modifies genetic risk for brain disorders. Analysis of post-mortem brain (n=304) revealed 8557 stress-interactive expression quantitative trait loci (eQTLs) that dysregulate expression of 915 eGenes in response to stress, and lie in stress-related transcription factor binding sites. Response to stress is robust across experimental paradigms: up to 50% of stress-interactive eGenes validate in glucocorticoid treated hiPSC-derived neurons (n=39 donors). Stress-interactive eGenes show brain region- and cell type-specificity, and, in post-mortem brain, implicate glial and endothelial mechanisms. Stress dysregulates long-term expression of disorder risk genes in a genotype-dependent manner; stress-interactive transcriptomic imputation uncovered 139 novel genes conferring brain disorder risk only in the context of traumatic stress. Molecular stress-encoding explains individualized responses to traumatic stress; incorporating trauma into genomic studies of brain disorders is likely to improve diagnosis, prognosis, and drug discovery.

Reduced LYNX1 expression in transcriptome of human iPSC-derived neural progenitors modeling fragile X syndrome.

Lack of FMR1 protein results in fragile X syndrome (FXS), which is the most common inherited intellectual disability syndrome and serves as an excellent model disease to study molecular mechanisms resulting in neuropsychiatric comorbidities. We compared the transcriptomes of human neural progenitors (NPCs) generated from patient-derived induced pluripotent stem cells (iPSCs) of three FXS and three control male donors. Altered expression of RAD51C, PPIL3, GUCY1A2, MYD88, TRAPPC4, LYNX1, and GTF2A1L in FXS NPCs suggested changes related to triplet repeat instability, RNA splicing, testes development, and pathways previously shown to be affected in FXS. LYNX1 is a cholinergic brake of tissue plasminogen activator (tPA)-dependent plasticity, and its reduced expression was consistent with augmented tPA-dependent radial glial process growth in NPCs derived from FXS iPSC lines. There was evidence of human iPSC line donor-dependent variation reflecting potentially phenotypic variation. NPCs derived from an FXS male with concomitant epilepsy expressed differently several epilepsy-related genes, including genes shown to cause the auditory epilepsy phenotype in the murine model of FXS. Functional enrichment analysis highlighted regulation of insulin-like growth factor pathway in NPCs modeling FXS with epilepsy. Our results demonstrated potential of human iPSCs in disease modeling for discovery and development of therapeutic interventions by showing early gene expression changes in FXS iPSC-derived NPCs consistent with the known pathophysiological changes in FXS and by revealing disturbed FXS progenitor growth linked to reduced expression of LYNX1, suggesting dysregulated cholinergic system.

Rescue of deficits by Brwd1 copy number restoration in the Ts65Dn mouse model of Down syndrome.

With an incidence of ~1 in 800 births, Down syndrome (DS) is the most common chromosomal condition linked to intellectual disability worldwide. While the genetic basis of DS has been identified as a triplication of chromosome 21 (HSA21), the genes encoded from HSA21 that directly contribute to cognitive deficits remain incompletely understood. Here, we found that the HSA21-encoded chromatin effector, BRWD1, was upregulated in neurons derived from iPS cells from an individual with Down syndrome and brain of trisomic mice. We showed that selective copy number restoration of Brwd1 in trisomic animals rescued deficits in hippocampal LTP, cognition and gene expression. We demonstrated that Brwd1 tightly binds the BAF chromatin remodeling complex, and that increased Brwd1 expression promotes BAF genomic mistargeting. Importantly, Brwd1 renormalization rescued aberrant BAF localization, along with associated changes in chromatin accessibility and gene expression. These findings establish BRWD1 as a key epigenomic mediator of normal neurodevelopment and an important contributor to DS-related phenotypes.